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The Kepler problem is a dynamical system that is well defined not only on the 
Euclidean plane but also on the sphere and on the hyperbolic plane. First, the theory 
Qh' of central potentials on spaces of constant curvature is studied. All the mathematical 

expressions are presented using the curvature k as a parameter, in such a way that 
they reduce to the appropriate property for the system on the sphere S 2 , or on the 
hyperbolic plane H 2 , when particularized for k > 0, or k < 0, respectively; in addition, 
the Euclidean case arises as the particular case k = 0. In the second part we study the 
main properties of the Kepler problem on spaces with curvature, we solve the equations 
and we obtain the explicit expressions of the orbits by using two different methods: first 
by direct integration and second by obtaining the K-dependent version of the Binet's 
equation. The final part of the article, that has a more geometric character, is devoted 
to the study of the theory of conies on spaces of constant curvature. 
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1 Introduction 



The hydrogen atom in a spherical geometry S* 3 was firstly studied by Schrodinger [1] in 
1940 and analyzed by Infeld [2] and Stevenson [3] at the next year; some five years later 
Infeld and Schild [4] obtained the spectrum of this system in an "open universe of constant 
negative curvature" (hyperbolic space H 3 ). Later on, Higgs [5] and Leemon [6] analyzed 
the characteristics of the two fundamental central potentials, Kepler problem and harmonic 
oscillator, on the iV-dimensional sphere and Kalnins et al, in a study of the dynamical 
0(4) symmetry of the hydrogen atom [7], considered the separation of variables of the 
Schrodinger equation with a Coulomb potentialn on the sphere S 3 using a classifications of 
coordinate systems separating the Laplace equation in the Riemannian spaces of constant 
curvature previously obtained by Olevskii [8] . Since then, a certain number of authors have 
studied this question from both the classical (Kepler problem) and the quantum (hydrogen 
atom) points of view [9] -[24]. The Kepler problem and the harmonic oscillator are the two 
most important superseparable systems and, because of this, their spherical and hyperbolic 
versions have been occasionally obtained, in some of the quoted references, as particular 
cases in the study of superintegrable systems on non- Euclidean configuration spaces. 

The Euclidean superintegrable systems in the plane were studied by Fris et al [25] and 
in the three-dimensional space by Evans [26]; later on different authors have considered 
this question [27]-[39] from different points of view. In Ref. [16] we studied the existence 
of the spherical and hyperbolic versions of these superintegrable systems taking the space 
curvature as a parameter; next the properties of the isotropic and nonisotropic harmonic 
oscillator and of those superintegrable systems related with the harmonic oscillator were 
analyzed [19, 20, 23, 24]. Now, our objective is the study of the Kepler problem on the 
sphere S 2 and the hyperbolic plane H 2 , both coming naturally from reduction from the 
'true' Kepler problem in a three-dimensional curved space S 3 or H 3 . 

In differential geometric terms, the Euclidean plane IE 2 is in fact a very particular 'limiting 
case' of the constant curvature spaces; accordingly, in dynamical terms, certain classical and 
well known potentials (Kepler problem, harmonic oscillator, and systems related with them) 
can also be considered as very particular 'limiting cases' of more general 'curved' systems. 
If the curvature k is introduced as a parameter then the question is the analysis of certain 
potentials in the space with constant curvature k (we may term these as ^-dependent po- 
tentials) with appropriate flat limit. As many different /t-dependent functions may have the 
same limit when K-»0we must require that certain fundamental properties of the Euclidean 
system continue to hold for the 'curved' system. By fundamental properties we mean those 
related with separability (superseparability) and integrability (superintegrability). The im- 
portant point is that the spherical and hyperbolic Kepler potential, to be studied in this 
article, can be characterized as a very specific /t-deformation of the well known Euclidean 
system, and conversely, the Euclidean Kepler potential arises as the particular k — case of 
the ^-dependent 'curved' system. 

The main purpose of this article is to study and solve the Kepler problem in the space of 
constant curvature k. As we will see, one of the advantages of this /t-dependent formalism is 
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that the study runs parallel to that of the Euclidean case. So the structure of this article is 
very similar to that of the corresponding chapter (or sections) devoted to central potentials 
and Kepler problem in some classical books of theoretical mechanics [40, 41]. In fact this is 
the idea of using the deformation as an approach: the curvature k will modify many things 
but will preserve the fundamental structure; this is so because the deep properties ensuring 
the possibility of exact solution for the Euclidean Kepler problem are preserved when a 
constant curvature is ascribed to the space. Of course, there are significant differences; one 
of them will be the geometric interpretation of the solutions that, although rather simple 
in the Euclidean plane, will pose some interesting geometric questions in spaces of constant 
curvature. 

In more detail, the plan of this article is as follows: In Sec. II we analyze the free geodesic 
motion on the spaces S 2 , IE 2 , and H 2 , as well as the existence of /t-dependent Noether 
symmetries. In Sec. Ill we study some general questions concerning the central potentials 
on spaces of constant curvature. Sec. IV, that can be considered as the central part of the 
article, is devoted to the study the ^-dependent Kepler problem on the spaces (S 2 , IE 2 , H 2 ). 
We have divided this section in three subsections: in the first part we consider the equivalent 
one-dimensional problem and we draw a rough classification of the orbits without requiring 
explicit solutions; in the second part we solve explicitly the problem and we obtain the 
expressions of the orbits; we arrive to the results by using two different methods: firstly by 
direct integration and secondly by obtaining the /t-dependent version of the Binet's equation. 
The third subsection is devoted to the analysis and classification of the ^-dependent orbits; 
first on the sphere S* 2 (k > 0), and then in the hyperbolic plane H 2 (k < 0). This third 
subsection poses the difficult problem of interpreting the ^-dependent equations as conies in 
curved spaces; of course all the obtained results reduce to well known Euclidean trajectories 
when we consider the particular value k — 0. Sec. V, that has a geometric character, is 
devoted to the study of (non- Euclidean) conies on the sphere S 2 and on the hyperbolic plane 
H 2 , which leads to the identification of the orbits of the Kepler problem with conies for 
any value of the curvature. Finally, in Sec. VI we discuss the results and make some final 
comments. 

2 Geodesic motion, Noether symmetries and Constants 
of Motion on (S 2 ,E 2 ,# 2 ) 

In the following we will make use of the same notation and techniques introduced for the 
oscillator in the above mentioned previous articles. That is, from now on all the mathematical 
expressions will depend of the curvature k as a parameter, in such a way that when assuming 
k > 0, k = 0, or k < 0, we will obtain the corresponding property particularized on 
the sphere, on the Euclidean plane, or on the hyperbolic plane. In order to present these 
expressions in a form which holds simultaneously for any value of k, we will use the following 
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"tagged" trigonometric functions 

C K (x) 



COS y/KX if At > 0, 

1 if k = 0, 

cosh^/— if « < 0, 



S«(a;) 



i 



sm v/ti 



-J= sinh 



if k > 0, 
if « = 0, 

■KI if At < 0, 



and the At-dependent tangent function T K (x) defined in the natural way, T K (x) = S K (x) / C K (x)- 
The fundamental properties of these curvature-dependent trigonometric functions are 



CI(x) + kS 2 k (x) 



and 



C K (2x) = 
S K (2x) 



Cl(x)-KSl(x) 



2S K {x) C K (x), 



dx Ck\X) , 

^C K (x) = -«S«(a:) 



(2) 



It is well known that the appropriate coordinates for the study of central potentials are 
polar (r, 0) coordinates (see Appendix); according to the intrinsic viewpoint, the coordinate 
r in any space of constant curvature At is the actual distance measured along the geodesies 
emanating from an arbitrarily chosen origin point O, while has the same meaning as in 
E 2 , the angle between this geodesic and a reference geodesic through O [42]. The range of 
r depends on At; specifically r G [0, ti/2^Jk\ for At > and r G [0, oo] for At < 0. The range of 
is the interval [0, 2tt]. So we start with the following At-dependent expression 

ds 2 = dr 2 + S 2 (r)d0 2 , 

that represents the differential line element on the spaces (S 2 , IE 2 , H 2 ) with constant curva- 
ture At. This metric reduces to 

ds{ = dr 2 + (sin 2 r) d(f) 2 , ds 2 = dr 2 + r 2 d(j) 2 , ds 2 _ l = dr 2 + (sinh 2 r) d(f) 2 , 

in the three particular cases of the unit sphere At = 1, Euclidean plane At = 0, and 'unit' 
Lobachewski plane At = — 1. 

The three At-dependent vector fields Yi, Y 2 , Yj, with coordinate expressions given by 



Yi(«) 

^(«) 

Yj 



(cos 

(sine; 

d_ 

00' 



dr 



q(r) 

S,c(r) 

q(r) 

S«(r) 



sm c 



cos 



00 ' 



are generators of three different one-parameter groups of diffeomorphisms preserving the 
metric ds 2 (isometries of the Riemannian manifold) [24]. In fact, the commutators of these 
vector fields are given by 



[Yi( K ) , Y 2 {n)\ =-kYj, [Yi(k) ,Yj] = ~ Y 2 (k) , \Y 2 (k) , Yj\ = Y^k) , 



4 



so that they close the Lie algebra of the group of isometries of the spherical (Euclidean, 
hyperbolic) space. Notice that only when k — (Euclidean plane) Y x and Y 2 will commute. 
Moreover, the Lagrangian for the geodesic (free) motion is given by the kinetic term arising 
from the Riemannian metric; that is 

L o (r,0,w r ,^;/t) =T(k) = (i) (w r 2 + S«(r)uj) , 

and is invariant under the actions of Yi(n), Y 2 (k), and Yj(k). 

A general standard Lagrangian (^-dependent kinetic term minus a potential) has the fol- 
lowing form 

L(r, 0, v r , k) = (-) ( v 2 + S 2 K (r)vl )-U(r, 0; «) , 
in such a way that for n = we recover a standard Euclidean system 

lim L = {\){ v 2 + r 2 v 2 )- V(r, 0) , V(r, 0) = U(r, 0; 0) . 

In some particular cases this Lagrangian system possesses the vector fields Y ± , Y 2 , or Yj, as 
exact Noether symmetries. If we denote by Yf, s = 1,2, J, the natural lift to the tangent 
bundle (phase space) of the vector field Y s and by 6l the Cartan semibasic one- form [43] 

dL dL 

6 L = —— dr + —— d(f> = v r dr + S«(r)% d<f> , 

OV r OVff, 

then the cases with exact Noether symmetries are the following: 

1. If the potential U is of the form U = U(z 2 ), with z 2 = S K (r) sin0, then 

Pi(«) =i(Yt(K))9 L = (cos0)t; r - (C«(r) S«(r) sin0)^ 
is a constant of motion. 

2. If the potential U is of the form U — U(z±), z\ — S K (f) cos0, then 

P 2 (k) =i(Y 2 t (K))9 L = (sin0)w r + (C«(r) S K (r) cos0) ^ 
is a constant of motion. 

3. If the potential C/ depends only on the coordinate r (central potential), then 

J(K) = t(Yj)9 L = sl(r)v^ 

is a constant of motion. 



We remark that the dependence on the curvature k is only explicit in the radial dependence 
of these functions; the angular dependence, contained in cos0 or sin0, is /t-independent. 
As an example, the vector field Yj is /t-independent but the integral of motion J(k) is in- 
dependent. These quantities P±, P 2 , J, are the /t-dependent versions of the two components 
of the linear momentum and the angular momentum. 
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3 Central potentials on (S 2 , E 2 , H 2 ) 



Making use of the three /t-dependent functions, P 1 , P 2 and J, the kinetic energy T(n) can 
be rewritten as follows 



T(k) = ( 1 -)(p?( K ) + P^) + kJ\k)), 

so that the total energy becomes 

E(k) = (\) ( P?(k) + + k J\k) ) + U(r, 0; k) , 

showing that, on spaces of (constant) non-zero curvature, the angular momentum has a 
contribution to the total energy of the system, proportional to the curvature k. 

A /t-dependent potential U whose expression in 'polar' coordinates (r, 0) has the structure 

G(0) 



U = F(r) + 



S*(r) ' 



is Hamilton- Jacobi separable [24]. It is, therefore, integrable and has the following two 
quadratic integrals of motion 



h{n) = P^) + P^) + 2F(r) + 2 ^, 
7 2 («) = J 2 (k)+2G(<P). 

Then, in this separable case, the total energy splits as a sum of two independent constants 
of motion 

= (i)/i(«) + (i)«/ 2 («). 

Of course, if G — then U = F(r) is a central potential and the function hi^) just reduces 
to I 2 = J 2 {k). 

The Lagrangian of a /t-dependent central potential is given by 

m = (l){v 2 r + sl(r)vl)-U(r-K), 

so that the dynamics is represented by the vector field 

d d . d . d 

XL = Vr d-r +V *dj + fr d^r +U d^ 

with the functions f r and given by 

/r = S„(r) C K (r>J-[/;, 
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The associated two equations are 

j t v r = S«(r) C K (r)vl -U' r , U' r = ^ , 
^J(k) = 0, J(«)=S 2 K (rK, 

where the time-derivative d/di can be interpreted, in geometric terms, as the Lie deriva- 
tive along Xl. The 0-equation just gives the conservation law of the re-dependent angular 
momentum (or the law of areas); concerning the radial equation it can be rewritten as follows 



d 

Vr = r = -dt 



1 J 2 



^S 2 K (r) 

so that, after multiplying by v r = r, it leads to 
that represents the conservation law of the energy 



o, 



^ = 0. 

dt 



Solving for r, we obtain 



dr 

di~ \ 



E - U(r; re) - 
that can be solved for dt and integrated 



,L J 2 



'2 y S 2 K (r) 



dr 



2[E-U(r;K)-(±)(jysUr))] 

so that in the particular re = case we recover the integral of the Euclidean case appearing 
in the books of theoretical mechanics 



t 



dr 



fi[E-V(r)-(i)(J*/r*)] 



Coming back to the general re ^ case, if we consider a change of variable from r to a new 
variable u K = C K (r)/ S«(r), by using (2) we get 



dr = 



du K 



ul + re 



and by considering the potential U — U(r; k) as a function of the new radial variable, 
U = U(u K ; k), the above integral for t becomes 



' J « + k) p [E p - U(u K ; k) - {\).Pul)] ' 
where E P denotes a constant of motion which which plays an important role: 

E p = E-{ 1 -)kJ 2 . 

(Notice the true energy for curvature n is E and not Ep] however Ep can be seen as a kind 
of /t-deformation of the Euclidean energy because E P — E for k — 0). This integral gives 
the value of t as a function of u K and consequently also of r; it can be inverted, at least 
formally, and we can obtain u K (or r) as a function of t and the two constants E P and J. 
The expression of as a function of t is just given by 

(f> = J J(u 2 K + K,)dt + (j) . 

We close this section by calling the attention to the factor (uf, + k) in the above integral for 
t that clearly resembles to the factor appearing in the elliptic integrals of the third kind. 



4 The Kepler problem on (S 2 , E 2 , H 2 ) 

The following spherical (hyperbolic) Lagrangian with curvature k, 
M«) = (l) {vr + Sl(r)v} ) - U K (r; «) , U K = - 



T«(r) 



represents the /t-dependent version of the Euclidean Kepler problem; the potential Uk re- 
duces to 

Ui(r) = --!?-, Uo (r) = V(r) = --, U^r) = - -\- , 
tanr r tanhr 

in the three particular cases of the unit sphere (k — 1), Euclidean plane (k = 0), and 'unit' 
Lobachewski plane (k = —1); the Euclidean function V(r) appears in this formalism as 
making separation between two different behaviours (see Fig. 1). The global sign has been 
chosen so that k > corresponds to an attractive potential. 

This potential is actually worthy of the name 'Kepler in curvature «' due to two reasons. 
First, this potential is the spherically symmetric potential satisfying Gauss law in a three- 
dimensional space of constant curvature k, where the area of a sphere of radius r is 4n S^r), 
whence the flux of the corresponding radial force field across a sphere of radius r, which 
clearly equals 4iiS 2 l (r)(dUK/dr) should be a constant independent of r; this condition leads 
directly to the potential Uk (As in the Euclidean case, the potential —(k/ T K (r)) does not 
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satisfy the Gauss law in two dimensions). The second reason is based on superintegrability: 
what singularizes Uk among other central potentials with the correct —k/r Euclidean limit 
is the property of being a superintegrable system for all the values of the curvature k [16]. 
In fact, Uk is endowed with the following two additional integrals of motion 

h(n) = P2(k)J(k) — k cosc/>, 
= Pi(k) J{k) + k sin , 

that represent the two-dimensional curvature versions of the Runge-Lenz constant of motion, 
whose existence is a consequence of the additional separability of Uk{t\ k) in two different 
systems of K-dependent 'parabolic' coordinates (this superseparability is not studied in this 
article). Of course only three of the four integrals are functionally independent; so we have 
several possibilities for the choice of a fundamental set of /t-dependent integrals of motion 
as, for example, {E(k), J(k), {E(k), J(k), I 4 (k)}, or {E(k), I 4 (k)}. 



4.1 The classification of orbits and the equivalent one— dimensional 
problem 

Let us start for an arbitrary central potential U(r). We have previously obtained 



so, if we introduce the 'effective one-dimensional equivalent ^-dependent potential W defined 

J 2 



as 

7-2 



W K (r) = U(r; k) + 



2S 2 (r)' 

where the term J 2 / (2 S«(V)) plays clearly the role of the 'centrifugal barrier potential', then 
the above property reduces to conservation of the energy for the fictitious /t-dependent one- 
dimensional problem arising from the potential W K 

(\)v 2 r + W K (r) = E 

so that we can use, also in this /t-dependent case, the method of the classification of the 
orbits by analyzing the behaviour of W K . 

Next we will study the main characteristics of the orbits for the Kepler potential in cur- 
vature k where W K is given by 

k T 2 



T«(r) 2S 2 (r) 
so it reduces to 

k J 2 

W 1 = - + —r-r- , 

tan r 2 sin r 
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k . J2 

W = 



r 2r 2 
W-x = - - — — + 



tanh r 2 sinh 2 r ' 

in the three particular one-dimensional problems associated to the unit sphere (k = 1), 
Euclidean plane (k = 0), and 'unit' Lobachewski plane (k = — 1). 

(1) Analysis of the potentials W\ and W c (c = ^/k > 0). 

The function W\ satisfies the following limits in the boundaries 

lim WAr) = +oo , lim WAr) = +oo , 

it cuts the r-axis at the points solutions of the equation sin(2ri j2 ) = J 2 /k, and it has a 
minimum at the point r s m given by r s m = tg _1 (J 2 /k). It represents, therefore, a potential well 
with barriers of infinite height at both extremes, r = and r = ir, and one single minimum 
placed inside the left half-interval (r s m < vr/2) with a value W s m = (l/2)( J 2 - k 2 / J 2 ) (Fig. 2). 
Thus, all the trajectories are bounded and all of them describe non-linear one-dimensional 
oscillations between the two turning points. 

In the general case of the function W C) the points r li2 are given by the two roots of 
sin(2cr 12 ) = c(J 2 /k) and the minimum is placed in the point r s m = (1/c) tg~ 1 (cJ 2 /A;) with a 
value given by = (l/2)(c 2 J 2 — k 2 / J 2 ). All these expressions are c-dependent and have 
the correct limits for Wq: 

J 2 J 2 k 2 
limri = — - , limr 2 = +oo , limr* = — , lim Wi, = . 

c^o 2k ' c-o ' c-o m k c-o m 2J 2 

(see Fig. 3). 

(2) Analysis of the potentials W-\ and W- c (c = v / — « > 0). 
The function W_i satisfies the two following limits 

lim W-i(r) = +oo , lim W-Ar) — — k , 

it cuts the r-axis in the point r\ unique solution of the equation sinh(2rx) = J 2 /k and, in 
the case that k and J satisfy the condition J 2 /k < 1, then it has a unique minimum in the 
point given by r h m = tanh _1 ( J 2 /k) with the value W^ = - (1/2) (J 2 + A; 2 /J 2 ). There 
exist therefore two possible situations: 



1. If J 2 /k < 1 the motion is bounded (periodic) for small energies W^<E < —k, and 
unbounded for higher energies E> — k. 

2. If J 2 j k>l the function W_\ will take the form of a potential barrier with infinite height 
at the the origin r = 0; only energies satisfying E> — k will be allowed and all the 
trajectories will be unbounded (scattering) open curves. 
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These two possible behaviours are represented in Figure 4. 
The general hyperbolic function W^ c satisfies 



lim W- C (r) = +00 , lim W_ c (r) = — ck , 



r->0 



it cuts the r-axis in the point r\ solution of sinh(2cri) = c(J 2 /k) and, if the condition 
c(J 2 /k) < 1 is satisfied, then it has a unique minimum in the point = (1/c) tanfT 1 c(J 2 /k) 
with a value given by = — (l/2)(c 2 J 2 + k 2 /J 2 ). It is clear that the smaller the value 
of c is, the easier of satisfying is the condition for the existence of a well, so that the behaviour 
(1) is becoming more and more dominant. Finally, the Euclidean limit is given by 

J2 J2 p 

lim r\ = — , lim = — , lim = 



2k ' "^o' m k ' c^6" m 2 J 2 
The convergence of W- c into IVo is represented in Figure 5. 

4.2 Determination of the orbits of the /^-dependent Kepler prob- 
lem 

4.2.1 Method I: Direct Integration 

We have previously obtained, making use of the conservation of the total energy E and the 
angular momentum J, two expressions for r and 0, which can be written as 

dr ^ ,S 2 K (r 



dt = . dt = ( 



{2 [E - U - (1/2) (J7S*(r))]' V J 

Eliminating t between both equations we have 

d , = Jdr 

Sl(r)y/2[E-U-(l/2) (J 2 /S 2 K (r))] ' 

that after the change of variable r — > u K with u K — \j T K (r), dr = — S 2 K (r) du K , becomes 
where-R(Mfe) denotes the following function 

„/ x 2E 2U , r, Q . 

= -j^-^-« + r 2 K ) 
2E P 2U 2 

" "j 2 ~J2 ~ U k- 
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All this is valid for a general potential. Next we particularize for the /t-Kepler problem, 
Uk{t\ k) — — kj T K (r), and then we obtain 

<f) = O du " , R(u k ) =a + f3u K + <yu 2 K , 



>R(u K ) 

with coefficients a, f3 and 7 given by 

2E P 2k 

In this particular case the integration is elementary and we arrive at 

_^ u K -(k/J 2 ) ^ r 2j2 W 

leading to 

n (rh\ = ( 

■J 2 



u K {4>) = (t^JI 1 + e K cos(0 - O )] , e K = y/l + z K . (3) 



This is the polar equation of the orbit in either space (S 2 ., IE 2 , if 2 ), for any value of the 
curvature k. It reduces to 



1 / k \ 

— = (^)[l + eicos(0-0 o )], 



tanr ^J 2 
1 / k \ 

~ = ( J2 ) t 1 + e ° C0S( ^ ~~ ^o)l ' 

^ = (-^)[l + e_ 1 cos(0-0 o )], 

in the three particular cases of the unit sphere (k — 1), Euclidean plane (/c = 0), and 'unit' 
Lobachewski plane (k — — 1) respectively. 

We recall that, as only three of the several integrals of motion can be functionally in- 
dependent, there must exist some relations between them. Concerning the couple h(K-), 
h(n), that, as stated above, represent the ^-dependent version of the Euclidean Runge-Lenz 
constants of motion, they are related with E P and J by 

II + I\ = 2E P J 2 + k 2 

so that we arrive to 

1 



6k = k ^ + 

that represents the natural extension to the curvature- dependent case of a well known prop- 
erty of the Euclidean case. 

Let us close this direct integration approach with three observations. First, it turns out 
that in any of the three manifolds (S 2 ., IE 2 , i/ 2 ), and for any value of k, this curve is a conic 
with a focus at the origin, where 'conic' has to be taken in a metric sense, relative to the 
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intrinsic metric in each space. This follows from the geometrical study to be presented in 
Section 5. Second, notice that the quantity e K is related to the 'partial energy' E P exactly 
as the Euclidean eccentricity of the conic is related to the total energy, yet only in the 
Euclidean case (k = 0) this function Ep coincides with E. In the general case (k 7^ 0) the 
characteristics of the orbit will be more easily associated to the values of the e K and/or E P 
than to the total energy E. Third, both the method and the results obtained show a close 
similarity with the Euclidean ones. In fact, the important point is that the classical and well 
known change of variable r—>u = 1/r admits as a generalization the /t-dependent change 
r^u K — 1/ T K (r) which affords a significant simplification for all values of k; conversely, this 
^-dependent change reduces to the Euclidean change r^u — 1/r for k — 0. 



4.2.2 Method II: Equation of Binet 



The definition of the angular momentum J determines a ^-dependent relation between the 
differentials of the time t and the angle 

Jdt = S 2 K (r)d(f). 

The corresponding relation between the derivatives with respect to t and is 

d / J \ d 
dt = ^sl(ry #' 

so that the second derivative with respect to t is given by 



d 2 



(; 



J ^ d 



(; 



J \ d 



dt 2 v S 2 (r) y rf0 LV S 2 (r) y rf0 . 

Introducing this notation in the radial equation, it becomes 



( 



J ^ d 



Sir 



J \ dr 



Sir 



\q3 r )J r 



Stir 



This equation can be simplified in two steps: firstly, the lefthand side can be rewritten by 
making use of 

d / C K (r) \ = / 1 \ dr 
dAs K (r))~ " V S 2 (r) J rf0 ' 

secondly, we introduce the change r — > u K in such a way that the potential U — U(r; k) be 
considered as a function of u K ; then we have 



U' 



dU 

du K . 



In this way we arrive at the differential equation of the orbit 

d 2 u K / 1 x 



d(p 2 
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that permits us to obtain as a function of u K for the given potential U (r; k) when considered 
as function of u K : 

# = /{«- (|)f-<}- (1/2) *«. 

Let us now particularize for the Kepler problem. In this case the potential Uk is given by 
Uk = — ku K , and the equation reduces to a linear equation with constant coefficients 

d 2 u K k 

which has the general solution 

k / k \ 

u K = A cos(0 - O ) + — = (_) [1 + e cos(0 - <j> )] , (4) 

where A (or e = A (J 2 /k)) and 0o are the two constants of integration. 

Remark that the differential equation of the orbit, usually known as the Binet's Equation, 
is preserved by the /t-deformation. That is, u deforms to u K but the the equation by itself 
remains invariant. 



4.3 Analysis of the orbits 

From the geometrical viewpoint, one of the integration constants in (3) or (4) can be made 
to disappear by simply choosing to be measured from the orbit position closest to the 
focus. Thus the orbit is 

T,(r) = — - - (5) 
1 + e cos (p 

which depends on two geometric parameters, the constants D and e. Comparing with (3) or 
(4), these two geometric parameters are related to the angular momentum J and the total 
energy E by means of the expressions 



^ J 2 2J 2 / J\ 

D = ^ e « = V 1 + l^- K T)' (6) 

remaining valid for any value of the curvature k and reducing to the known expressions for 
the Euclidean = 0. Thus, as advanced before, the relation between the geometric 

constants D, e, and the physical constants J, E — (1/2) kJ 2 = E P) is the same for all values 
of curvature; here we see again the quantity Ep to play an role: it can be considered as a 
kind of 'translational' part of the energy, not taking into account the contribution to the 
energy of the angular momentum. 

The nature of the orbit depends on the values of these constants. We discuss separately 
the cases with positive and negative k. 
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(i) Spherical space, k > 0. The polar equation of the orbits is: 




(7^)[ 1 + eK COS(0-0 O )]- 



Tfc(r) tan^r 



The orbit is always a closed curve and for any value of e K this curve is an spherical ellipse 
with a focus at the origin. For a fixed J the minimal value of the total energy corresponds 
to orbits with e K = 0, which are circles with radius r = J 2 /k and total energy E cir = 
(l/2)(—k 2 /J 2 + kJ 2 ) (which may be either negative or positive). The possible energies for 
given J fill the interval [E C \ T) oo] and correspond to the parameter e K in the interval [0, oo] 
and to Ep in the interval [— (l/2)(k 2 / J 2 ), oo]. In the following description the focus will 
be conventionally placed at the North pole, and upper and lower half-spheres refer to the 
North and South hemispheres. This orbit has always a vertex closest to the focus given by 
the unique solution of the equation 1/ T re (r) = {k/ J 2 ){1 + e K ), and depending on the value 
of e K we may however distinguish three possible behaviours for the second vertex: 

• When e K < 1 the orbit is completely contained in the upper half-sphere centered at the 
focus, and the values of r always remain less than the length 7t/(2 v /k) of a quadrant on 
the sphere of curvature k. In some respects this reminds the case of Euclidean ellipses. 

• If e K — 1, the ellipse has the second vertex at r = it/ {2^/k) (reducing to r = n/2 when 
k — 1), and is a curve touching tangentially the equator associated to the focus at the 
origin. When k — > 0, the limit of the distance r = n/(2y/K) between the origin and 
the equator on a sphere S 2 is r = oo and thus this curve is analogous to a Euclidean 
parabola, touching the spatial infinity at a point. 

• When e K > 1, the orbit has the second vertex at r > 7t/(2 v /k); it is a big spherical 
ellipse crossing the equator of the focus and entering into the lower half-sphere. This 
is somehow analogous to an Euclidean hyperbola. 

(ii) Hyperbolic space k < 0. Let us consider the general case of arbitrary negative curvature. 
The orbit equation is now 



and the relations (6) also apply but now with k < 0. For a fixed value of J, the diagram of 
the effective equivalent potential reveals two essentially different situations for the behaviour 
of W(r), corresponding to values of J smaller or larger than the particular value Jqo of J: 



When J < Jqo (J 2 /k < 1/\/—k), the equivalent potential is qualitatively as its Euclidean 
counterpart; it has a minimum at r = r cir determined by T K ( r cir) = J 2 /k, with a negative 





k 
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value for W(r cir ) = E cir = —k 2 /2J 2 + kJ 2 /2 and when r increases from r cir to r — > oo, the 
effective potential tends from below to a negative constant value —ky/—n. When J = J^ 
then the minimum occurs at r = oo. When J > there is no longer a minimum in the 
equivalent potential, which is a monotone function of r, tending when r — > oo to the value 
—ky/—K from above. This discussion alone suffices to establish two landmark values for the 
energies (at a fixed value of J): 

• E = Ea V = (1/2) (—k 2 /J 2 + k.J 2 ), with J 2 /k < corresponding to the circular 
orbits with E Pcir = -(1/2) {k 2 / J 2 ), 

• E = E^ = —ky/^K, with Epoo = —ky/^K — (l/2)/cJ£, = —(l/2)kyT I K for orbits 
where the particle barely reaches spatial infinity with velocity 0. 

To go further in any analysis about the exact nature of the orbits is easier starting from 
some geometric knowledge of the properties of conies, to be sketched in the next section. 
However, at this point we can provide a first classification into closed and open orbits, and 
we may discuss its main characteristics. 

When the energy is in the interval [E C1[ , E^), the motion is bounded and periodic, while 
for E in the interval (E^, oo], the motion is not periodic and the orbit goes to the infinity; 
notice that E = is contained within this second interval, and thus is not the separating 
value between bounded and unbounded orbits. The angular momentum provides a second 
constant of motion, which may have values in the two subintervals J 2 / k e [T K (r c i r ), l/y/—n\ 
or J 2 /k G [l/y/—K, oo]; it is pertinent to discuss these two possibilities successively. 

Let us first discuss the subcase J < J^, where the equivalent potential has a minimum. 

1. If the parameter e K = then the solution is clearly a circular orbit (circle in the 
hyperbolic plane) of radius r cir such that T K (r cir ) = J 2 /k which happens precisely at 
the minimum for the equivalent potential. The total energy for a circular orbit is given 
by E cir = - {\/2){k 2 /J 2 - k.J 2 ) that corresponds to E PaT = -{\/2){k 2 / J 2 ). 

This is, the angular momentum of any circular orbit should remain lower than the limit 
value Joo which is approached for the circular orbit with r — > oo. In the Euclidean 
case the value of J is not bounded for circular orbits (Joo tends to infinity for k — > 0). 

2. If e K is greater than but in the interval < e K < 1 — \J — k (J 2 /k) then the orbit 
will be a closed curve not reaching the spatial infinity. As we will check in Sec. 5, this 
curve is indeed an hyperbolic ellipse, with a focus at the potential origin. The interval 
of allowed values for the total energy E turns out to be 




or equivalently, if we consider the "translational" part of the energy: 
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3. If e K is in the remaining interval, 1 — \J— k (J 2 /k) < e K , corresponding to the values 
—\J—k k < E < oo of the energy E, then the the motion is non periodic and the orbits 
are unbounded open curves. The border value for e between these closed periodic and 
open not periodic types of orbits is the value e K = 1 — \J —k (J 2 /k), which should be 
thought of as corresponding to the 'last ellipse', or horoellipse. It is not immediately 
clear whether the orbits with e above this value are parabolas or hyperbolas in the 
hyperbolic plane. However, there is natural to expect some special status for the 
orbits with E = 0, e K = \Jl — /t(J 4 //c 2 ) > 1, which we might expect to be parabolas 
in the hyperbolic plane. This would imply that at least all orbits in a finite interval of 
values of the e K parameter, whose lower bound is 1 — \J —k (J 2 /k) and which contains 

e K = \J\ — k (J 4 /k 2 ) are parabolas; notice this interval should reduce to the single 
value e K = 1 in the Euclidean limit k — > 0, so the existence of many different parabolic 
orbits does not spoil the known fact that there is only a parabolic orbit in E 2 . All 
these conjectural properties are indeed true, as we will confirm in Sect. 5. 

We can summarize these three points as follows: the different types of trajectories in H 2 , 
k < as a function of the parameter e K are the following 

hyperbolic circle = e K 

hyperbolic ellipses < e K < 1 — \J — n (J 2 /k) 

hyperbolic horoellipse 1 — \J — k (J 2 /k) = e K 

hyperbolic open conies 1 — \/—k (J 2 /k) < e K < oo . 

Alternatively, the landmark values can also be given for the classification as a function of 
the total Energy E: 

hyperbolic circle E = — (1/2) (k 2 / J 2 - kJ 2 ) 

hyperbolic ellipses - (l/2)(/c 2 /J 2 - kJ 2 ) < E < -ky/^K 

hyperbolic horoellipse E = — k \J — k 

hyperbolic open conies — k \/—k < E < oo 

Thus the behaviour of the hyperbolic dynamics shows some important differences with re- 
spect to the Euclidean case. We postpone until the next section the establishment of the 
particular value of e K , or E, separating the open conies into parabolic and hyperbolic regimes, 
and we simply state that e K = \Jl — k (L 4 /k 2 ), corresponding to E — 0, gives an orbit which 
is a parabola in H 2 . The value of e K , or the corresponding energy E, separating from 
parabolic to hyperbolic trajectories is a /t-dependent positive value to be determined; this 
will be done in the next section. When the parameter e K is greater than the upper value for 
parabolic orbits, then we will have hyperbolas in H 2 as orbits. 

The remaining cases happen when J > J^. For these values of J the equivalent potential 
is a decreasing function of r with no minimum at all, and hence the possible motions should 
have energies above E^; these would be always non-periodic motions and open orbits. We 
will classify them into parabolas and hyperbolas in the next section. 
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5 Conies on spaces of constant curvature 



In this section we give a geometric description of conies in the three k > 0, k = 0, k < 0, 
constant curvature spaces S 2 , E 2 , H 2 , emphasizing those aspects relevant in relation with 
Kepler motion in these spaces. To stay within a reasonable extension limits we restrict 
to a mainly 'statement of facts' presentation, which should either confirm the conjectures 
advanced in the dynamical part, or serve as a geometric foundation for them. The important 
thing is the perfect matching between the dynamical approach in Sec. 4 and the results here. 

As customary in this context, lines will mean the geodesies of the constant curvature space. 
The metric intrinsic geometric definition of conies, that can be applied to any 2-d space of 
constant curvature k involves focal elements, i.e., either oriented points or cooriented lines. 
In any such space, and by definition: 

An ellipse/hyperbola will be the set of points with a constant sum/difference 2a of distances 
7*1, r 2 , to two fixed points F 1 ,F 2 , called foci and separated a distance 2f. 

A parabola will be the set of points with a constant sum/ difference 2a of distances r\, f 2 , to 
a fixed point F\, called focus, and to a fixed cooriented line / 2 , called focal line; the oriented 
distance 2p> between iq and / 2 plays here the role of focal separation. 

An ultraellipse/ultrahyperbola will be the set of points with a constant sum/ difference 2a 
of oriented distances fi,f 2 , to two fixed intersecting lines /i,/ 2 , separated by an angle 2F 
and called focal lines. 

In the generic case of constant curvature k ^ these three pairs of curves, each pair sharing 
the same focal elements, are the generic conies; Euclidean plane is not generic among the 
family of constant curvature spaces, thus some Euclidean properties of conies are very special 
and do not provide a good starting viewpoint to discuss the k ^ properties. 

Further to generic conies, particular distinguished conies appear for instance when the focal 
separation vanishes (/ = 0, ip — 0, F — 0) and limiting conies when some focal elements go 
to infinity (if possible at all); both particular and limiting conies can be obtained as suitable 
limits from the generic ones. 

We first list some particular conies. When the focal separation vanishes / = 0, the ellipse is 
a circle and the hyperbola is a pair of intersecting lines. No useful intuition about parabolas 
can be drawn from the non-generic Euclidean case; a parabola with vanishing focal separation 
ip — is only a particular instance of parabola, to be named here equiparabola, but when k ^ 
there are other parabolas, with any positive or negative <p ^ 0, which are not equiparabolas. 
The main point is that general parabolas are not to be considered here as limiting conies, 
as usually done in E 2 . An ultraellipse with focal angle F = is a ultracircle (a equidistant 
curve, the set of points equidistant from a fixed line), while an ultrahyperbola with F = 
is a pair of lines. 

The three generic types of conies, as well as these particular conies with zero focal separa- 
tion do exist in any space of constant curvature k. Further to that, when the curvature is 
negative, there are new types of limiting conies, where either the foci or the focal lines may 
go to infinity. 
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5.1 Conies in geodesic polar coordinates in a space of curvature k, 

Now let us describe conies in the three essentially different spaces S 2 , E 2 , H 2 ,. 

In the sphere S 2 , with k > the situation is rather simple. First, as there are neither 
points nor lines at the infinity, there are no limiting cases in the sense they will appear in H 2 . 
Second, as there is a I- 1 correspondence (the usual polarity) between oriented points and 
cooriented lines (North Pole versus Equator), with a constant distance tc/2^/~k separating 
them, it follows that any focal element (either point or line) implies the existence of a 
polar focal element (line or point); thus any conic of one of the three given type pairs may 
be considered as a conic of any other prescribed type pair. And further, as the distance 
between two antipodal points is constant tt/^/k,, by changing a focus F 2 to its antipodal 
point F 2 , (which according to the definition is again a focus of the conic), an hyperbola 
with focus Fi,F 2 , will be as well an ellipse with focus F 1: F 2 . A similar possibility exists 
for focal lines which may be changed to their antipodal lines, with coorientation changed. 
Thus all spherical conies can be seen as spherical ellipses, including the two particular zero 
focal separation cases, i.e., circles (spherical parallel circles) and pairs of lines (a pair of 
intersecting large circles in the sphere). For the purposes of getting a common view we 
may still think in three generic types, keeping in mind that all three coincide and there is 
complete freedom in understanding a given spherical conic as either an ellipse, or a parabola 
or a hyperbola (just as a geodesic circle can be seen at the same time as a geodesic circle, 
with center at O or as an equidistant curve to the polar of O); more on this will be said 
later. 

In the Euclidean plane E 2 ultraellipses and ultrahyperbolas (with any focal angle F) are 
just pairs of parallel lines, in directions parallel to the two bisectors of the focal angle; in the 
standard algebraic classification of Euclidean conies these appear as degenerate conies. But 
when compared to S 2 or H 2 , another more important degeneracy happens in E 2 : two families 
of parabolas with the same focus and any two parallel focal lines will coincide (because a pair 
of parallel lines are equidistant in E 2 ). This means that the parabolic focal separation between 
the focus and the focal line can be chosen arbitrarily for a fixed Euclidean parabola, which 
only determines its focus and the direction of the focal line. As a consequence parabolas 
are no longer generic conies in E 2 . While parabolas have always some intermediate status 
between ellipses and hyperbolas, in the Euclidean case parabolas appears exclusively as 
limiting cases between ellipses and hyperbolas, and not as the full fledged species of conies 
they are in the generic k ^ case. 

The hyperbolic case H 2 with negative constant curvature is richer. The limiting conies 
obtained from a ellipse/hyperbola through a point P when a focus stays fixed and the other 
goes to infinity are called horoellipse/horohyperbola. The limiting conies obtained from a 
ultraellipse/ultrahyperbola through a point P when a focal line stays fixed and the other 
goes to infinity are called horoultraellipse/horoultrahyperbola. These four types of conies are 
also limiting forms of parabolas, with either focus fixed and focal line going to infinity, or 
focal line fixed and focus going to infinity. In E 2 all these limiting conies collapse precisely 
to either Euclidean parabolas or pairs of lines. 
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There are more limiting conies in H 2 : the horocycle, obtained either from a circle through 
P when its center goes to infinity, or from an ultracircle through P when the baseline goes 
to infinity, and the osculating parabola, the limit of a parabola when the focus go to infinity 
along the fixed focal line. 

In order to make contact with the results in the dynamical part, let us draw our attention 
precisely on the conies with a proper focus. From now on we completely disregard ultrael- 
lipses and ultrahyperbolas, which have two focal lines, and osculating parabolas, which have 
not a proper focus. We consider exclusively ellipses/hyperbolas and parabolas, as well as 
their limiting cases — horoellipses and horohyperbolas — and particular cases — circles, lines 
and equiparabolas corresponding to vanishing focal separations — . All these conies have a 
proper focus, which we may place at the origin point, as well as another focal element, which 
may be another focus or a focal line. In either case the conies have as a symmetry axis the 
line joining the two foci or going through the focus and orthogonal to the focal line. To 
match with the familiar Euclidean expressions, we introduce polar coordinates in S 2 , E 2 , 
H 2 , taking the origin at the fixed focus and the symmetry axis of the conic as the half-line 
= 0. With this choice, the common equation of the whole family of conies described above 
may be written in the form: 

T„(r) = — - - (7) 
1 + e cos 

with D and e non-negative. This can be derived by trigonometric considerations from the 
definition of the conic in the space of constant curvature k using the relations given in [44] 
and will be discussed elsewhere. This relation justifies the claim made in Sec. 5: the orbits 
of the Kepler problem in either S 2 , E 2 , H 2 , are conies with a focus at the origin of the 
potential. The dependence on the polar angle is exactly the same as in the Euclidean 
case and this may suggest to consider e as the non-Euclidean analogue of the eccentricity 
of the conic, though when k ^ the link between the values of e and the type of the conic 
is not so direct as in E 2 . We should mention that there are at least two other quantities 
keeping different properties of the Euclidean eccentricity, and no single quantity keeps all 
properties; in this paper we will only be concerned with e, and for brevity we shall refer to 
it as eccentricity. 

For an ellipse/hyperbola with focal distance 2f and 2a as the sum/difference of distances to 
the focus, and for a parabola with focal separation 2tp and 2a as sum/difference of distances 
to focus and focal line, the eccentricities turn out to be 

S*(2/) C K (2yQ 
6ell/hyP " S«(2a)' 6par " C K (2a) 

reducing in the Euclidean k = case to the well known f /a for ellipses/hyperbolas and 1 
for parabolas. 

In the polar coordinate system (r, 0), when k > or k = 0, the range of values of T K (r) is 
the whole real line (completed with oo), but when k < 0, the values of T«(r) are confined to 
the interval [0, or to [—l/y/—K, 1/V~ K \ if negative values for r are allowed along 

the opposite semi-axis according to the usual practice. 
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Now, for the conic (7) the minimum value for r happens when = 0. This means that D 
and e are independent constants when k > but must fulfil the inequality < D/{1 + e) < 
1 / \J — k when k < 0. 

Notice that in (7) the periastron of the orbit is placed on the semi-axis = at r = r per , 
and thus r per is related to D and e by 

D 

T nypcr) 



1 + e 

which for any value of k has always a unique root for r per . Next, let us look for the intersection 
of the conic with the line orthogonal to the conic symmetry axis through the focus, which 
corresponds to = tt/2. The distance between this intersection point and the focus is 
traditionally called semilatus rectum of the conic (in the Euclidean case) and will be denoted 
by p. For k — all Euclidean conies (except the limiting double straight line with p = oo) 
intersect this semilatus line. For any k ^ 0, p would satisfy 

T«(p) = D, 

but there is a difference between the cases with k > or k < 0. While for k > 0, any 
D will determine a unique p, when k < this equation will define a real semilatus rectum 
only when D < 1/ \/—k, the semilatus rectum will be formally infinite when D = 1/ yj—n 
and will not exist at all when D > This means that in the case of the hyperbolic 

plane H% the family of conies we are considering includes conies intersecting the 'semilatus' 
line at a proper point, only at infinity or not intersecting at all. We have arrived to the 
following situation: all conies with given fixed focus and symmetry axis intersect this line in 
exactly one does not intersect in E 2 and an infinite number of them do not intersect in 

HI 

Thus, when k > the equation of the complete family of conies we are considering is 

T„(r) = - ^ {P) - (8) 
1 + e cos 

but in the negative curvature case k < 0, there will be some conies — precisely those with 
D > 1 / \f--K — not described under (8). To cater for these cases it will prove useful to 
introduce another real distance p, complementary to the ideal semilatus rectum, related to 
Dby 

7 \ 7^ = D 

(-«) T K (p) 

and in terms of this choice, when k < the equation of the complete family of conies we are 
considering is given by one of the two mutually exclusive possibilities 

Tk ^ = l + eSs0' Tk ^ = (- K )T«(p)(l + ecos0) ' (9) 
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where the ranges are < p < oo and oo > p > respectively. The conic 

T ^ r ) = 7=71 1 ^ (10) 

\J—k (1 + ecos^J 

is the common limit p — > oo and p — > oo of (9). Notice the two expressions (9) can be used 
as well when k > but then each of the two alternatives covers actually all cases, and are 
then redundant. Only the first possibility in (9) has a sensible Euclidean limit because the 
k — > limit of the p family gives a straight line at the infinity of the Euclidean plane; thus 
the family p as a set of conies different from the p family is specific to the hyperbolic plane. 

Now the only remaining problem is to link the values of the two parameters D and e to the 
type the conic belongs. When is such a conic an ellipse, a parabola or a hyperbola according 
to their definitions in the space of curvature k! When is it a particular conic, either a circle, 
a equiparabola or a line? And finally, — only for the hyperbolic plane H 2 — , when is it a 
limiting horoellipse or horohyperbola? 

In the Euclidean plane the answer is well known and easy: the conic type depends only on 
e, and not on p. This is not so when k ^ 0. 



5.2 Conies on the sphere S% 

A given conic in S 2 will have a unique and well defined parameter e in (8). But labelling 
it as an ellipse, or parabola or hyperbola, requires to choose a particular set of two focal 
elements, and a ± sign to decide between the sum or difference in the definitions. Unlike 
on H 2 where these choices are unique, for the sphere they may be made in several ways. 
By means of suitable choices, any spherical conic with any < e < oo may be considered 
as an ellipse, as a parabola or as a hyperbola. This is an unavoidable consequence of the 
definitions, which are clearly natural ones for non-zero curvature. 

If we want to classify conies in S 2 into three disjoint species, this would require adopting 
an additional convention (somehow bypassing the definition), which may be chosen in several 
reasonable ways, but which nevertheless remains as a convention. 

A possibility would be to class spherical ellipses with < e < 1 as 'ellipses', those with 
e = 1 as 'parabolas', and those with e > 1 as 'hyperbolas'. This is the convention adopted by 
Higgs [5] who formulated it according to the property of not crossing, touching or crossing the 
equator, which can be easily shown to be equivalent to the requirement e<l,e = l,ore>l, 
respectively. Another different convention follows from the fact that any ellipse/hiperbola 
with focus Fi,F 2 (and focal separation 2f) can also be considered as a hyperbola/ellipse 
with focus F 1 , F 2 (and focal separation 2f = 2(n/2y/K, — /)). The sum of the two focal half- 
separations / and / is 7r/2y / K, this is, the length of a quadrant on the sphere. Thus, we may 
consider as 'ellipses' exclusively those spherical ellipses with a value of the focal separation 2f 
less than a quadrant, that is, 2f < 7r/2y / K, and as 'hyperbolas' exclusively those spherical 
ellipses with a value of 2f greater than a quadrant, that is, 2f > n/2y/R. Within this 
convention, and if we want artificially enforce non-redundancy, the name 'parabolas' will 
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only be left for the spherical ellipses with precisely 2f = 7r/2y / K, which in our previous 
nomenclature were called equiparabola (as 2f = ix /2y/R between two focus means that one 
focus will be incident with the the polar of the other focus, which is the focal line). In this 
view 'ellipses' are completely contained in the half-sphere with center at the ellipse center, 
while 'hyperbolas' will go through the boundary of the half-sphere centered at the hyperbola 
center. 



5.3 Conies on the hyperbolic plane H\ 

The case of the hyperbolic plane H\ is different and there is no conventionality in it. In the 
following paragraphs we are implicitly assuming k < 0. 

When is the conic (9) an ellipse? In addition to the actual vertex at the point of closest 
approach to the focus, placed on the = axis, ellipses will have another vertex, placed on 
the = 7i semi-axis with a positive value for T K (V a po) which anyhow should belong to the 
interval [T K (V pcri ), \/y/—K,\. The lower bound in this interval corresponds to circular orbits 
while the upper places the apoastron at the infinity, and the ellipse will go to an horoellipse. 
It is easy to conclude that the positivity condition requires e to be in the interval [0, 1] and 
then T K (r apo ) G [T K (r pcri ), l/yj—n\ only happens within the first alternative in (9), provided 
that the eccentricities lie in the interval 

< e e u < 1 - v 73 ^ T«(p) • 

Notice this interval depends on p and for any p ^ is strictly smaller than the Euclidean 
one, which is recovered of course, when k — > 0. The lower bound for e corresponds to circular 
orbits, while the upper bound is the limiting value corresponding to horoellipses 

e cir = 0, e ho roell = 1 - T K (p) . 

When will the conic (9) be an hyperbola! In this case, in addition to the vertex at the 
periastron, the conic will have another vertex point (on the other hyperbola branch, thus 
not actually on the physical orbit) on the = semi-axis, which will appear on the equation 
as a negative value for T K (r apo ) on the = n semi-axis. Hence, the values of T K (r) for 
= 7r must belong to the interval [—1/^—k, — T K (r pcr )], with the lower value corresponding 
to horohyperbolas, and the upper value to lines. By using the same strategy, it is easy to 
conclude that the negativity condition requires e G [1, oo], and then T K (r apo ) may lie in the 
required interval in either of the two alternatives (9), with e parameters respectively in the 
intervals 

1 + T K (p) < e hyp , 1 + ^_ - — < e hyp , 

Notice both conditions lead to a common limit when p, p — > oo, namely, 2 < eh ypC rboias for 
the conies (10). Here particular and limiting conies correspond to: 

1 

Chorohyp 1 V ^ Tk\P) i ^hoj-ohyp 1 1 ~ j Clincs . 

y— K \-k\P) 
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Once we have characterized the intervals of eccentricity for ellipses and hyperbolas, the 
remaining gap correspond to parabolas. Thus, within the first family: 

1 - T«(p) < e par < 1 + T«(p) , 

and within the second 

1 \ <e 

The particular parabolas with zero focal separation, this is, the equiparabolas, corresponds 
to the values: 

e e quipar = 1/ Ck(p) = \j 1 + KT«(p), e cquipar = Cfc(p) = 1/^1 + KT«(p) • 

and are different from the parabolas with e = 1. 

This way we get the following result: conies (9) in the Hyperbolic plane H 2 , with a fixed 
value of p may be either ellipses, parabolas or hyperbolas, according to the value of their 
eccentricity as follows 

circle e C j r = 

ellipses < e e ii < 1 — \J — k T k (p) 

horoellipse e hor oeii = 1 - v^-k T K (p) 

parabolas 1 - T K (p) < e par < 1 + y/^n T K (p) 

horohyperbolas e h oroh yP = 1 + T K (p) 

hyperbolas 1 + T K (p) < e hyp < oo 

and conies (9) with a fixed p may be only parabolas or hyperbolas: 

The most remarkable property is the existence, for a fixed value of the semilatus rectum p, 
of a full interval of values for eccentricity corresponding to parabolas, including always two 
special values, e — 1/ Chip) an d e = 1; in the Euclidean limit k — > this interval collapses 
to the single value e par = 1, and there is a single Euclidean parabola for a given p. 

Now let us consider the complete family of conies corresponding to trajectories with a 
given periastron distance, say r per ; these would correspond to the trajectories of a particle 
launched orthogonally to the radial direction from a point at distance r per from the potential 
origin and with a given velocity. The corresponding pattern for the Kepler problem in the 
Euclidean case is well known: a circle, a family of ellipses of increasing eccentricity < e < 1, 
a separating parabola, a family of hyperbolas of increasing eccentricity 1 < e < oo and finally 
a limiting straight line corresponding to infinite velocity. How about the analogue of this 
pattern in the non-zero curvature case? 

First we require all conies will have a periastron at r pcr . Within the two families the 
relation between p or p with r per is given by 

T K (r pcr ) = ^ T K (r pcr ) = ^—^—^ 
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and, using these expression in the classification stated above, we get the landmark values of 
p or p corresponding to the different types of conic: 



which in the Euclidean limit, with k — and To(x) = x, reduces to the well known: 



Figures 6A and 6B represent a set of Kepler orbits in the hyperbolic plane H 2 . 

6 Final comments and outlook 

We have solved the Kepler problem on the three spaces of constant curvature. As we have 
stated in the introduction, one of the fundamental characteristic of this approach is the 
use of the curvature k as a parameter. In this way, all the ^-dependent properties that we 
have obtained reduce to the appropriate property for the system on the sphere S 2 , or the 
on hyperbolic plane H 2 , when particularized for k > 0, or k < 0, respectively; in addition, 
the Euclidean case arises as the very particular (but important) case k — 0. So, we can 
summarize this situation pointing out two important facts: 

• The Kepler problem is not a specific or special characteristic of the Euclidean space 
but it is well defined in all the three spaces of constant curvature. 

• There are not three different Kepler problems but only one that is defined, at the same 
time, in three different manifolds. 

Of course, since the three manifolds are geometrically different, many dynamical properties 
show differences according to the characteristics of the manifolds; nevertheless, the important 
point is that there is only one theory that is simultaneously valid for the three manifolds 
and for any value of the curvature. We illustrate this situation recalling the following two 
important results: 

• The equation of Binet remains true in the three cases. 

• The orbits are conies in the three cases. 




r pcr < Pcii < 2r per = p par = 2r per < p hyp , when k = . 



For the hyperbolic p family, which has no Euclidean limit, we have: 
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It is well known that the classical Euclidean Kepler problem is one of the systems endowed 
with more interesting properties. Therefore, it would be interesting to study, all of them, 
in the general non-Euclidean k ^ case; as an example, we have made use, in Sec. 4, 
of two constants of motion, I^k) and I^k), that represent the K-dependent version of the 
Runge-Lenz vector. The Runge-Lenz vector plays a very important role in the Euclidean 
case, so it is natural to suppose that the same situation will be true for the case of the general 
'curved' system. In fact, we think that, as the Euclidean system is just a very particular 
case of a much more general system, all the Euclidean characteristics, that are obtained 
and discussed in the books of theoretical mechanics, must admit a ^-dependent deformed 
version appropriate for the general "curved" system. We think that these questions are 
open problems that must be investigated. We also note that in all this paper we have only 
dealt with the two-dimensional case since we assumed that it arises from a reduction of the 
three-dimensional Kepler problem on S 3 and if 3 ; this dimensional reduction is worthy of a 
further study. 

Finally, the analysis of the orbits has led us, in a natural way, to the theory of conies on 
spaces of constant curvature. Notice that, although Sees. 3 and 4 were mainly concerned with 
dynamical questions, Sec. 5 was written emphasizing its geometrical character; in fact, it can 
be considered by itself and independently of the other previous sections. It is clear that the 
theory of conies on spaces of constant curvature is a geometrical matter of great importance 
deserving a deeper study that we hope to present elsewhere. On the one hand, there exist 
some general points and basic properties that, although perfectly known and clearly stated 
in the Euclidean plane, remain to be studied in the case of spaces of constant curvature; 
on the other hand, we have already obtained some particular points as, for example, the 
existence several different parabolas in the Hyperbolic plane, or the possibility of some 
different alternative ways of defining the eccentricity for a conic in a k ^ space, which are 
really noteworthy. These geometrical questions are also open problems to be investigated. 

Appendix: Polar geodesic coordinates 

A two-dimensional manifold M can be described by using different coordinate systems. If we 
consider it as an imbedded submanifold of 1R 3 , then the points of M can be characterized by 
the three external coordinates, as (x,y,z) or (r,(f>,0), plus an additional constraint relation. 
Nevertheless, in differential geometric terms, a more appropriate approach is to develop the 
study by using two-dimensional systems of coordinates adapted to M. 

On any general two-dimensional Riemannian space, not necessarily of constant curvature, 
there are two distinguished types of local coordinate systems: "geodesic parallel" and "geo- 
desic polar" coordinates. They reduce to the familiar cartesian and polar coordinates on 
the Euclidean plane (see Refs. [16] and [42]) and both are based on a origin point O and a 
oriented geodesic g± through O. 

For any point P in some suitable neighbourhood of O, there is a unique geodesic g joining 
P with O. The (geodesic) polar coordinates (r, 0) of P, relative to the origin O and the 
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positive geodesic ray of g±, are the (positive) distance r between P and O measured along 
g, and the angle between g and the positive ray g±, measured around O (Fig. 7). These 
coordinates are defined in a neighbourhood of O not extending beyond the cut locus of O; 
polar coordinates are singular at O, and <p is discontinuous on the positive ray of g 1 . 

In the case of M being a space of constant curvature k, the expression for the differential 
element of distance ds 2 is given by 

ds^dr' + sl^d^, 

so that we get ds 2 = dr 2 + r 2 d(p 2 for the particular k — Euclidean case. 
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ure Captions 



Figure 1. Plot of the Kepler Potential as a function of r, for the unit sphere 
k — 1 (upper curve), Euclidean plane k = (dash line), and 'unit' Lobachewski plane 
k — — 1 (lower curve). The three functions are singular at r = but the Euclidean 
function Uq — V appears in this formalism as making a separation between two different 
behaviours. In fact U = V is the only Potential that vanish at long distances. 

Figure 2. Plot of W c as a function of r, for k — 1 with k and J given by (k — 1, J = 
!)• 

Figure 3. Plot of W c as a function of r, for (/c = 1, J = 1) and three different values 
of the curvature c = : c — 1.2 (upper curve), c = 0.8 (middle curve), and c = 0.4 
(lower curve). 

Figure 4. Plot of W- c as a function of r, for c = 1 (/c = — 1). The upper curve 
corresponds to (k — 1, J = 2) and lower curve corresponds to (A; = 4, J = 1). 

Figure 5. Plot of M^_ c as a function of r, for (k — 4, J = 1) and three different 
values of the curvature c = ^/ — k : c = 1.5 (lower curve), c = 1 (middle curve), and 
c = 0.5 (upper curve). 

Figure 6A and 6B. A set of Kepler orbits with a fixed periastron distance, depicted 
in the conformal Poincare disk model of hyperbolic plane, for 'small' r pcr (Figure 6A) 
and 'large' r per (Figure 6B). The potential center is at the origin, which is a focus 
of the conies. Thick lines represent particular and limiting conies: circle, horoellipse, 
horohyperbola and straight line. A suitable selection (with the other focal elements 
chosen as to make the diagram clear) of ellipses, parabolas and hyperbolas are rep- 
resented as thin lines in each of the three ranges determined by the previous conies. 
For ellipses and hiperbolas the other focus (not marked) is on the horizontal line; for 
parabolas the focal line is perpendicular to this horizontal line. The semilatus rectum 
of the conic lies on the vertical straight line through the origin; notice in the figure 
(A) only some hyperbolas do not intersect this line, while in (B) conies which do not 
intersect the semilatus line include all hyperbolas as well as many parabolas; these are 
the two generic behaviours, as explained in the text. 

Figure 7. Polar geodesic coordinates in a two-dimensional Riemannian space (M, g) 
(Figure 7A) and in the sphere S 2 (Figure 7B). 
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U1,V,U_1 




Figure 1. Plot of the Kepler Potential as a function of r, for the unit sphere k = 1 
(upper curve), Euclidean plane k — (dash line), and 'unit' Lobachewski plane k — — 1 
(lower curve). The three functions are singular at r = but the Euclidean function Uq = V 
appears in this formalism as making a separation between two different behaviours. In fact 
Uq = V is the only Potential that vanish at long distances. 
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Figure 2. Plot of W c as a function of r, for k — 1 with k and J given by (k = 1, J = 1). 




Figure 3. Plot of W c as a function of r, for (/c = 1, J = 1) and three different values of 
the curvature c = y/n : c = 1.2 (upper curve), c = 0.8 (middle curve), and c = 0.4 (lower 
curve). 
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Figure 4. Plot of W_ c as a function of r, for c = 1 (k = — 1). The upper curve 
corresponds to (k — 1, J = 2) and lower curve corresponds to (k = 4, J = 1). 




Figure 5. Plot of jy_ c as a function of r, for (/c = 4, J = 1) and three different values 
of the curvature c = \J — k : c = 1.5 (lower curve), c = 1 (middle curve), and c = 0.5 (upper 
curve). 
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Figure 6 A and 6B. A set of Kepler orbits with a fixed periastron distance, depicted 
in the conformal Poincare disk model of hyperbolic plane, for 'small' r per (Figure 6A) and 
'large' r per (Figure 6B). The potential center is at the origin, which is a focus of the conies. 
Thick lines represent particular and limiting conies: circle, horoellipse, horohyperbola and 
straight line. A suitable selection (with the other focal elements chosen as to make the 
diagram clear) of ellipses, parabolas and hyperbolas are represented as thin lines in each of 
the three ranges determined by the previous conies. For ellipses and hiperbolas the other 
focus (not marked) is on the horizontal line; for parabolas the focal line is perpendicular 
to this horizontal line. The semilatus rectum of the conic lies on the vertical straight line 
through the origin; notice in the figure (A) only some hyperbolas do not intersect this line, 
while in (B) conies which do not intersect the semilatus line include all hyperbolas as well 
as many parabolas; these are the two generic behaviours, as explained in the text. 
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Figure 6 A and 6B. A set of Kepler orbits with a fixed periastron distance, depicted 
in the conformal Poincare disk model of hyperbolic plane, for 'small' r per (Figure 6A) and 
'large' r per (Figure 6B). The potential center is at the origin, which is a focus of the conies. 
Thick lines represent particular and limiting conies: circle, horoellipse, horohyperbola and 
straight line. A suitable selection (with the other focal elements chosen as to make the 
diagram clear) of ellipses, parabolas and hyperbolas are represented as thin lines in each of 
the three ranges determined by the previous conies. For ellipses and hiperbolas the other 
focus (not marked) is on the horizontal line; for parabolas the focal line is perpendicular 
to this horizontal line. The semilatus rectum of the conic lies on the vertical straight line 
through the origin; notice in the figure (A) only some hyperbolas do not intersect this line, 
while in (B) conies which do not intersect the semilatus line include all hyperbolas as well 
as many parabolas; these are the two generic behaviours, as explained in the text. 
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Figure 7 A. Polar geodesic coordinates in a two-dimensional Riemannian space (M,g). 
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